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ABSTRACT 



Aims. We derive the spectra of the cosmic X-ray background (CXB) and of the Galactic ridge X-ray emission (GRXE) in the 
~20-200keV range from the data of the IBIS instrument aboard the INTEGRAL satellite obtained during the four dedicated Earth- 
occultation observations in early 2006. 

Methods. We analyze the modulation of the IBIS/ISGRI detector counts induced by the passage of the Earth through the field of 
view of the instrument. Unlike previous studies, we do not fix the spectral shape of the various contributions, but model instead their 
spatial distribution and derive for each of them the expected modulation of the detector counts. The spectra of the diffuse emission 
components are obtained by fitting the normalizations of the model lightcurves to the observed modulation in different energy bins. 
Because of degeneracy, we guide the fits with a realistic choice of the input parameters and a constraint for spectral smoothness. 
Results. The obtained CXB spectrum is consistent with the historic HEAO-1 results and falls slightly below the spectrum derived with 
SwiftfBAT. A 10 % higher normalization of the CXB cannot be completely excluded, but it would imply an unrealistically high albedo 
of the Earth. The derived spectrum of the GRXE confirms the presence of a minimum around 80 keV with improved statistics and 
yields an estimate of ~0.6 M© for the average mass of white dwarfs in the Galaxy. The analysis also provides updated normalizations 
for the spectra of the Earth's albedo and the cosmic-ray induced atmospheric emission. 

Conclusions. This study demonstrates the potential of INTEGRAL Earth-occultation observations to derive the hard X-ray spectra of 
three fundamental components: the CXB, the GRXE and the Earth emission. Further observations would be extremely valuable to 
confirm our results with improved statistics. 

Key words. Earth ~ Galaxy: disk - Galaxies: active - diffuse radiation - X-rays: diffuse background - X-rays: general 



1. Introduction 

Although the cosmic X-ray background ( CXB) was discovere d 
before the cosmic microwave background (iGiacconi et al 
it is known in much less detail and its spectral shape and nor- 
malization are still subjects of debate. This diffuse emission 
is thought to be mainly due to unresolved active galactic nu- 
clei (AGN) extending to cosmological distances with a contri- 
bution from Type la supernovae in the low-energy gamma-rays 
JZdziai'ski 1996). Evidence for the AGN origin of the CXB at en- 
ergies below 10 keV comes from various X-ray mirror telescopes 
- in particular Chandra and XMM-Newton - that were able to 
reso lve up to 80% of the diffuse emission in to point sources 
(e.g. lBrandt & HasingedlloollGilli et al.ll2007h . The amount of 
resolved sources decreases rapidly with energy though so that 
only 2.5 % of the diffuse background is resolved by the deep- 
est surve y yet in the 20-60 k eV range, at the peak of the CXB 
emission dPaltani et alj|2008l) . The characterization of the actual 
spectral shape and normalization around this emission bump is 
crucial to estimate the fraction of heavily absorbed Compton- 
thick AGN tho ught to contribute significantly in t his hard X-ray 
spectral range dUeda et al.ll2003l; iGilli et al.ll2007l; ISazonov et al] 



l2b08HTreister et alJl2009h 

The High Energy Astronomical Observatory 1 {HEAO-1) is 
hitherto the only satellite which had a dedicated mechanism to 



disentangle the CXB from the instrumental background. By us- 
ing a movable Csl crystal with a thickness of 5 cm to cover 
part of the field of view (FoV), the HEAO-1 observations of the 
mid-1970s are still the most accurate and reliable measurements 
of the CXB spectral shape in the hard X-ray s dMarshall et al.l 
1 9801: iKinzer et al J 1 1 9971: ICiruber et al.1 fT999li . Without such a 
masking mechanism in recent space missions, a practical way 
to study this diffuse hard X-ray emission is to use the Earth as a 
screen occulting part of the background sky. For pointing satel- 
lites in low orbits around the Earth, our planet often crosses part 
of the field of view during normal operations. Such events can 
be analy zed in detail to evalu ate the CXB spectrum. This was 
do ne by Frontera et alj (l2007h for the BeppoSAX mission and 
bv'Aiello etalJ (l2008h for the data of the Burst Alert Telescope 
(BAT) aboard the Swift spacecraft. 

With its eccentri c three-days orbit, the INTEGRAL satellite 
dWinkler et al.l2003 h is close to the Earth only during the perigee 
passage when the instruments are not operating because of ex- 
cessive background in the radiation belts. In order to study the X- 
ray background, a series of four dedicated observations were per- 
formed in January and February 2006. The Earth was allowed to 
pass through the FoV of the instruments shortly after radiation- 
belt exit while the s pacecraft was aimed t o point towards a fixed 
position in the skv. iChurazov et alJd2007l) described these obser- 
vations by all four instruments aboard INTEGRAL in great detail 
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Fig. 1. Model IBIS/ISGRI images of the sky 
with (right) and without (left) instrumental vi- 
gnetting effects (see Fig.O. They show the ge- 
ometry and relative intensity (on a logarithmic 
scale from dark and blue to red and bright ar- 
eas) of the various components during the first 
EO at IJD = 2216.04 as derived for channel 3 
(~27keV). The bluish circle shows emission 
of the Earth occulting the diffuse sky back- 
ground (purple), the Galactic ridge (red strip), 
and the point sources (bright dots). The instru- 
mental background is ignored for clarity. The 
images extend over the partially coded FoV 
(28.8° X 29.2°). Point sources are convolved 
with a circular Gaussian typical for the instru- 
mental resolution ((T = 0.1°). 



and studied how the passage of the Earth modulates the detector 
counts by occulting part of the CXB. 

The difficulty of using the Earth to shield the CXB comes 
from the fact that the Earth is not dark in the hard X-rays. The 
emission from the Earth in the 20-200 keV range consists of 
two major contributions: the reflection of the CXB by the at- 
mosphere and its Compton emission under the bombardment 
by cosmic rays (CR). Disentangling the CXB occul t ation f rom 
the Earth emission is challenging. IChurazov et al ] (I2007h as- 
sumed the spectral shape of the CXB and of the two Earth emis- 
sion components and fitted their normalizations to the observed 
amplitude of the Earth modulation i n the data. The st udies of 
the BeppoSAX/PDS mea surements ( j Frontera et al.ll2007[) and the 
SwiftfBACT observations (lAiello et al] 20081) alsorelied on a pri- 
ori assumptions on the spectral shape of the CXB and the Earth 
emission. 

We present here a completely different approach for the 
analysis of the s ame I NTEGRAL observations as were used by 
IChurazo v et alJ (|2007), but only focusing o n the data of the 
IBIS/ISGRI instrument (Ubertini et alJl20()3l) . Instead of fixing 
the spectral shapes of the CXB and the Earth emission compo- 
nents, we aim to derive them based on a detailed modeling of the 
spatial distribution of these components and of all instrumental 
effects. If the Earth surface brightness differs significantly from 
the uniform CXB occultation, it is possible to disentangle these 
components based on the recorded modulation when the Earth is 
crossing the FoV. This has the potential to simultaneously derive 
the spectral shape of the CXB and of the Earth emission from 
the observations. 

Another difficulty of the analysis is the presence of the 
Galactic plane in the border of the wide FoV of IBIS (see Fig.[T]). 
An empty extragalactic field would have been ideal to study the 
CXB, but this was not possible due to various scheduling con- 
straints. This complication is however an opportunity to study 
in a ddition the diffiise Galactic ridge X-ray emission (GRXE) 



ay . 

(e.g. Re vnivtsev et alJ l2006': Kri vonos et al teoOllBouchet et alJ 
[2008), which is also occulted by the passage of the Earth. 

The observational material and the analysis method are de- 
scribed in Sects. |2] and [3] respectively. We present the obtained 
spectra in Sect. |4] and discuss them in comparison to previous 
results in Sect.|5] before concluding in Sect.|6l Unless otherwise 
stated, the quoted errors are 1-cr uncertainties, i.e. at the 68 % 
confidence level (CL). 



2. Data 

The data used here are the four Earth-occultation observations 
(EOs) conducted by INTEGRAL in January and February 2006 
at the start of satellite revolutions number 401, 404, 405 and 406. 
We ref er the reader to the det ailed description of these observa- 
tions in lChurazov et al ] (l2007l) . We focus our analysis on the data 
of the IBIS/ISGRI gamma-ray imager that is best suited to study 
the emission in the ~20-200keV range. 

Our work is based on the analysis of the modulation in full 
detector lightcurves induced by the passage of the Earth through 
the FoV. These lightcurves are obtained with the latest version 
of the ii .light executable that will be included in a forth- 
coming release of the Off-line Scientific Analysis (OSA) soft- 
ware p ackage provided by the I NTEGRAL Science Data Centre 
(ISDC, 'Courvoisi er et al.ll2003h . They are corrected for instru- 
mental dead time and the effect of dead and noisy pixels, which 
amount to typically 5 % of all ISGRI pixels. For each of the four 
similar observations, we extracted detector lightcurves with a 
time binning of 300 s in a series of 16 energy bins (see Table|2]i, 
carefully chosen to isolate instrumental emission features, in 
particular the broad lines at 26 and 3 1 ke V from CdTe and the 
naiTow lines at 59 keV from W, 7 5-77 keV from Pb and 82- 
84keV from Bi (iTerrier et aUlIOOl . 

The detector lightcurves originally expressed in units 
of countcm^^s"' were multiplied by 0.5(128 x 0.4cm)^ = 
1310.72 cm^, which is the area of the detector assumed by the 
standard ISGRI ancillary response file (ARE), describing the en- 
ergy dependence of the effective area of the instrument. The fac- 
tor 0.5 refers to the fraction of open coded mask elements, and 
0.4 cm is the size of the 128 x 128 ISGRI detector elements. 
Apart from this change of unit to have full detector lightcurves, 
the only other manipulation of the data was a cleaning of the 
lightcurves. This was done by removing points with uncertain- 
ties of more than twice the average uncertainty and by iteratively 
removing a few isolated outstanding points lying more than 3 cr 
away from the smoothed lightcurve with a smoothing window of 
30 minutes (i.e. 6 time bins). 

In order to subtract from the lightcurves the contribution of 
point sources in the FoV, we performed an image analysis sep- 
arately for the four Earth observations. This was possible since 
the drift of the satellite was moderate despite the absence of star 
trackers during the pointings. The image analysis was done in 
a standard way with the default background maps of OSA 7.0. 
We searched for all sources previously detected by ISGRI in- 
cluding the new sourc e IGR J 17062-6 143 already reported by 
IChurazov et al.l (l2007h . We then selected all sources that were 
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Table 1. List of point sources detected in at least one of the four 
EOs. 



Source name 


RA Dec 

( deg )" 


EOl 

r 10 


E02 

"erg 


E03 
cm 


E04 

1 -,6 
) 


IGR J14471-6319 


1 O 1 

zzl.ol 


Tn 

-63. zy 


- 


- 


- 


1.9 


IGR J 145 15-5542 


TTi on 
ZZZ. o9 


^o 

-55.00 


18.4 


- 


- 


- 


IGRJ14532-6356 


IIS. Si 


iiQ nQ 
—63.93 


- 


- 


- 


2.2 


IGR J 15094-6649 


ZZ/.J6 


oo 

— 60.8z 


4.1 


- 


- 


4.4 


PSRB 1509-5 8 


OTO AO 

ZZS.4B 


cn 1 A 
—59.14 


7.2 


7.5 


- 


5.6 


IGRJ15359-5750 


zji.y / 


C7 OQ 

-5 /.83 


2.9 


- 


- 


1.4 


H 1538-522 


ZJ3.00 


CO Qn 


14.2 


22.8 


- 


25.0 


4U 1543-624 


ZJ0.9B 


— OZ.5 / 


- 


4.5 


4.4 


- 


IGR 116167^957 


244.16 


^ n no 


- 


2.8 


- 


7.1 


IGR J 16207-5 129 


lA'z in 
z43.iy 


^ 1 '^n 


2.6 


6.4 


- 


- 


SWIFT J1626.6-5 156 


z4o.OJ 


c 1 n/1 
-51.94 


12.4 


11.1 


10.1 


8.3 


IGR J16283-4838 


lAn r\A 
z4 / .U4 


—46.00 


8.2 


- 


- 


- 


IGR J 163 18-4848 


247.95 


-48.82 


37.7 


- 


11.6 


7.1 


IGR J 16320-4751 


248.01 


-47.87 


8.1 


- 


- 


12.8 


4U 1626-67 


248.07 


-67.46 


13.0 


12.3 


8.9 


8.7 


IGR J 16377-6423 


249.57 


-64.35 




1.1 






IGR J 16393-4643 


249.77 


-46.70 




11.9 




1.6 


H 1636-536 


250.23 


-53.75 


7.1 


2.3 






IGR J 17008-6425 


255.20 


-64.43 








0.9 


XTE J1701-462 


255.24 


-46.19 




26.8 






GX 339^ 


255.71 


-48.79 


9.5 


28.7 


40.9 


44.3 


IGR J 17062-6 143 


256.57 


-61.71 








4.5 


NGC 6300 


259.25 


-62.82 


2.5 


3.1 


3.4 


2.8 


ESO 103-35 


279.58 


-65.43 






8.3 





Notes. source catalog position in right ascension (RA) and 
declination (Dec). **' derived model fluxes in the 20-50 keV band for 
each EO. 



detected with a significance of more than 2 o" in the 22-60 keV 
band. We chose this low significance threshold to minimize the 
contribution of the known point sources to the GRXE and the 
CXB. Sources with even less significance are more likely to be 
spurious and their global contribution will mostly cancel out 
with fake negative sources. We tested both a simple power- 
law and a bremsstrahlung model to fit the data. We found that 
for most sources the bremsstrahlung model gives a better phe- 
nomenological description of the data than a powerlaw because 
many sources have a convex spectral shape in this energy range. 
The fluxes derived for the brightest (> 3 cr) sources in our sam- 
ple are listed in Table [T] 

3. Method 

The detector lightcurves described above were modulated by the 
passage of the Earth through the FoV of IBIS. Our approach was 
to model these observations in detail to derive the expected mod- 
ulation of the detector counts for each emission component on 
the sky. This resulted in a series of model lightcurves in dif- 
ferent energy bins for each emission component and each of 
the four EOs. We then fitted the normalizations of these model 
lightcurves to the observed detector counts to derive the actual 
contribution of the diffuse emission components. 

This method requires the knowledge of the spacecraft posi- 
tion and attitude with respect to the Earth and to the background 
sky at any time, a description of the spatial distribution on the 
sky of the various emission components and also an accurate de- 
scription of the IBIS/ISGRI instrumental characteristics. These 
aspects are described in the three subsections below. The gener- 
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Table 2. Numerical values of the obtained spectra shown in 
Fig.|6] 



E 


AE 


^sky 




^ ear 








( keV )" 




( keV cm ^ 


s"' sr"' 


)* 




21.14 


1.44 


jy.j J 


2.53 


7.05 


1 

1.00 


1 c riQ 


CM 


24.01 


1.44 


4Z. Jy 


1.92 


8.63 


1 An 

1 .4 / 


1 Q Q/1 


KJ.JJ 


27.36 


1.91 


40.04 


1.81 


11.23 


1 ^ 1 


1 Q 1/^ 
IV. ZD 


4.0(S 


31.19 


1.91 


4Z.VU 


1.44 


13.79 


1 1 

I. II 


1 A AO 
10. 02 


Z. / J 


35.02 


1.91 


A 1 nc 


2.45 


17.02 


1 An 


1/1 1^ 
14. 


4.Z0 


38.85 


1.91 




2.55 


20.31 


1 cs 
1.00 


1110 
IZ. IZ 


Z.UJ 


43.16 


2.39 


JJ.4Z 


2.42 


23.87 


z.4d 


V.oU 


U.o4 


47.95 


2.39 


38.67 


1.95 


26.82 


2.16 


9.70 


1.35 


52.73 


2.39 


36.24 


2.42 


29.30 


2.00 


8.33 


1.40 


58.00 


2.87 


30.88 


3.41 


30.85 


2.14 


7.72 


1.37 


63.74 


2.87 


28.86 


2.81 


31.18 


2.48 


7.17 


1.78 


71.88 


5.27 


28.92 


3.34 


32.47 


3.11 


6.46 


1.31 


83.37 


6.22 


30.01 


6.35 


33.37 


3.23 


6.17 


0.80 


94.86 


5.27 


27.12 


4.89 


34.01 


2.71 


6.64 


0.84 


112.57 


12.45 


24.82 


5.94 


35.04 


3.90 


6.68 


0.99 


162.35 


37.34 


30.65 


9.55 


37.85 


7.46 


7.93 


2.86 



Notes. Central energy, E, and half-width, AE, of the bins. **' Fluxes, 
F, and statistical uncertainties, AF, for the sky background (sky), the 
Earth (ear) and the Galaxy (gal). On average over the region 
320° < / < 340° and|fc| < 5°. 



ation of the model lightcurves is described in Sect. 13.41 whereas 
the actual spectral fitting procedure is the subject of Sect. 13.51 

3.1. Satellite position and attitude 

To construct the images of the sky corresponding to each of the 
four Earth observations as illustrated in Fig. [T] we needed to 
know the exact attitude of the satellite and its distance to the 
Earth at any time. This information can be extracted from aux- 
iliary files provided by the mission operation centre (MOC) in 
Darmstadt. It was used to compute the position of the Earth cen- 
ter, the position of the geographic and magnetic poles, and the 
apparent radius of the Earth as a function of time, all expressed 
in degrees, using the IBIS/ISGRI instrument coordinates (F, Z). 
The spacecraft was close enough to the Earth at the beginning 
of the observation for the planet's sphericity to slightly affect 
its apparent radius. This was properly taken into account, as 
well as the ~ 100 km of obscuring atmosphere in the hard X- 
rays mentioned bv lChurazov et al.r(l2007 ). The magnetic pole in 
the Northern hemisphere is set to its 2005 position of 82.7° N, 
1 14.4° W. 

3.2. Spatial distribution of components 

Although we were only interested in the temporal modulation 
of counts on the full detector area, we needed a sufficiently pre- 
cise description of the spatial distribution of the emission com- 
ponents. We chose to define all of them by analytical functions 
that are described below and are illustrated in Fig. [T] 

The simplest component is the CXB that we assumed to be 
uniform on the sky. Although there is evide nce that the CXB has 
some large-scale inten sity variations (e.g. iBoughn et al. I I2OOI 
iRevnivtsev et al.l2008l) . they are small in amplitude ( i 2 %) and 
it would be very difficult to evaluate and account for a possible 
non-uniformity so close to the Galactic bulge. 

The GRXE was modeled with two perpendicular Lorentzian 
functions aligned with the Galactic coordinates. The full-width 



4 



M. Tiirler et al.: INTEGRAL hard X-ray spectra of the cosmic X-ray background and Galactic ridge emission 



at half maximum (F WHM) of the Lorentzians are of 21° 
(iKrivonos et al.ll2007l Fig. 7) and 1.2° ('Revni vtsev et al.ll2006i 
Fig. 5) respectively along the Galactic longitude, and lati- 
tude, b. The Lorentzian's maximum are at the Galactic center 
wi th a slight latitude displ acement of b - -0.15° as measure d 
bv lRevnivtsev et al.l (l2006h . As shown bv lKrivonos et al.l (l2007l) . 
this distribution matches well the COBEfDWSE map at 4.9yum, 
which was used by Bouchetet al. (2008) as a template for the 
GRXE below 120keV. 

We took special care to define the spatial distribution of the 
Earth's emission. There are two distinct compon ents to be taken 
into account: the CXB reflection by the Earth (Churaz ov et al. 
2008| ) and the CR-induced atmospheric emission (Sa zonov et al. 
2007h . The authors of these studies performed Monte-Carlo sim- 
ulations to derive both the spectrum and the surface brightness 
of the Earth emission. We used the latter results as a precise de- 
ter mination of the expecte d image of the Earth at hard X-rays. 

IChurazov et all (l2008l) found that the X-ray albedo of the 
Earth is limb-darkened at lower energies and limb-brightened at 
higher energies. At energies below ~100keV - where this com- 
ponent dominates the Earth emission - the emission is found 
to be limb-darkened, but slightly less than for a sphere emitting 
black-body radiation. Such an object would have a linear depen- 
dence of the flux with s cos 9, where is the zenith angle, i.e. 
the angle between the line-of-sight and the normal to the surface. 
For the Earth albedo below ~100keV they instead found an an- 
gular dependence of the reflected flux that can be approximated 
by Fiji) oc (1 -0.5 fi). We used this equation to define the Earth 
albedo component. 

We modeled the CR-induced emission of the Earth's atmo- 
sphere according to ISazon ov e t al. (2007, Eq. (7)). By setting 
the solar modulation potential to (f>- 0.5 - corresponding to the 
solar minimum during the EOs of 2006 - we can simplify this 
equation as: 

c^^l(\+^i)[\ + {R,J3.2fy^\ (1) 

where /i is as defined above and /?cut is the geomagnetic cut- 
off rigidity. In the dipole approximation of the Earth's magnetic 
field, the latter depends ma inly on the geomagnetic latitude 
as /^cut - 14.5 cos'* GV dSmart & Sh ea 2005 ). The resulting 
atmospheric emission of the Earth is a combination of relatively 
strong limb-darkening from the yU-dependence in Eq. ([T]i with en- 
hanced emission at the magnetic poles from the /Im-dependence. 

We note that we took into account for both Earth emission 
components the distortion of the surface brightness related to 
the fact that only a portion of the Earth's hemisphere can be seen 
when the spacecraft is relatively close to the planet. 

3.3. Instrumental characteristics 

As we wanted to fit real detector lightcurves with model 
lightcurves we needed to take into account the instrumental char- 
acteristics of the telescope in the modeling. This does not include 
detector responses, but all effects attenuating the incoming pho- 
ton field on its path from outside the telescope until reaching 
the detector plane. We identified five effects that affected the de- 
tector illumination depending on the direction of the incoming 
radiation and sometimes on its energy. The most obvious effect 
is the attenuation due to the coded-mask elements which block 
out about half of the incoming radiation. The second effect is the 
non-uniform exposure map which is caused by a partial illumi- 
nation of the detector for a source outside of the fully coded FoV. 
Another vignetting effect is induced by the Nomex honeycomb 




Fig. 2. Surfaces representing the five IBIS/ISGRI vignetting ef- 
fects that affect the incoming radiation until it reaches the de- 
tector plane. The effects are those corresponding to channel 12 
(~72keV) and are: (a) the IBIS tube transparency to off-axis ra- 
diation, (b) the energy independent exposure map, (c) the effec- 
tive coded-mask transparency, (d) the transmission of the Nomex 
structure supporting the mask, and (e) the absorption of the 
ISGRJ spider beams. The total vignetting effect is the product 
of the effects (b) to (e) with the addition of effect (a), which ex- 
tends well outside the partially coded FoV. 

structure that supports the coded mask. The two last effects are 
due to the IBIS/ISGRI spider beams separating the eight detec- 
tor modules and to the lead shielding of the IBIS telescope tube. 
The aluminum spider results in opacity at the lowest energies, 
and the lead shielding of the tube becomes transparent at the 
highest energies. 

The five effects mentioned above are included in the standard 
IBIS/ISGRI software for image reconstruction and spectral ex- 
traction, but as we worked directly with the detector lightcurves, 
we needed to account for these effects independently. Their mod- 
eling as used in this work is described below and is illustrated in 
Fig.|2] while the overall vignetting effect is shown in Fig.[T] 

The coded-mask transparency is ideally of 0.5 since there are 
as many elements open as closed. This is, however, only true at 
the center of the FoV. For off-axis sources there is an additional 
attenuation due to the thickness of 16 mm of the mask elements, 
which project a wider shadow on the detector for increasing off- 
axis angles. As the elements are made of tungsten - a strongly 
absorbing material - it is fair to assume the elements to be com- 
pletely opaque in the energy range considered here. As we were 
only interested in the net effect over the full detector plane and in 
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a simple analytical description we approximated the mask pat- 
tern as a giant chessboard of 46 equally-sized square elements 
on a side of 1064 mm. We then properly computed the additional 
shadow from radiation that crossed the border of the mask ele- 
ments and did not fall onto the shadow of other elements. 

The exposure map of IBIS/ISGRI is basically very simple 
with a value of 1 in the fully coded FoV and a linear decrease 
to zero in the partially coded FoV, except in the corners of the 
image where the decrease is quadratic. When we modeled this, 
we properly took into account the disposition of the eight mod- 
ules of the IBIS/ISGRI detector and the two-pixel wide space 
between them. 

The Nomex structure supporting the coded mask of IBIS is 
absorbing part of the incoming photons. This is corrected for in 
the OSA software by off-axis efficiency maps depending on en- 
ergy. In 2006, at the time of the EOs, these maps were still an 
approximation with only a dependence on the off-axis angle. We 
used the new maps introduced in the OSA 6.0 release that do in- 
clude an additional azimuthal dependence due to the alignment 
of the walls of the hexagonal tubes that the honeycomb structure 
is made of and also a correction for the tubes pointing ~ 0.5° 
away from the center of the FoV; a misalignment likely due to 
on-ground manipulations of the spacecraft. We note that the at- 
tenuation by the cosine of the off-axis angle is included in these 
maps. 

The eight IBIS/ISGRI detector modules are separated by an 
aluminum structure called the ISGRI spider. It is made of one 
beam along the Z axis and three perpendicular beams. According 
to the IBIS experiment interface document part B (EID-B) the 
beams have a trapezoidal section with a height of 48 mm a base 
of 9.5 mm and a wall angle of 4° resulting in an upper width of 
2.8 mm. The wall angle ensures that the spider is not casting a 
shadow in the fully coded FoV. However, in the partially coded 
FoV, the spider can mask up to two rows of ISGRI pixels. We 
modeled this effect in detail for each energy bin based on the 
corresponding attenuation length of Al. We found that for some 
specific off-axis directions, the ISGRI spider can result in an at- 
tenuation of the radiation on the detector plane of up to ~ 5 % in 
the lower energy bins. 

The last instrumental effect we considered is the IBIS tele- 
scope tube transparency. At higher energies, the tube becomes 
transparent to radiation from outside the fully coded FoV, giving 
an additional contribution to the detector lightcurves. The tube 
is made of two vertical walls perpendicular to the Y axis and 
of two walls inclined with an angle of 3.47° transverse to the Z 
axis. The walls are shielded with glued lead foils. The thickness 
of the Pb sheets for each wall is reported in Table 3.2.8.1 of the 
EID-B (issue 7.0). We used the values for the four upper sheets 
of the wall that are relevant for the calculation of the tube trans- 
parency for an off-axis angle up to about 45°. The calculation 
was done carefully, avoiding radiation further blocked by other 
parts of the ISGRI collimating system: the 1 mm thick W shield 
of the ISGRI hopper and the 1 .2 mm thick W-strips of the side 
shielding of the mask. The tube transparency is neglectable at the 
lower energies, but reaches ~ 1 % in channel 12 (see Fig.|2]i near 
the Pb attenuation length edge at about 75-80 keV and ~ 6 % in 
the last energy bin. Although this might seem unrelevant, the cu- 
mulated effect on the detector plane from a wide region outside 
the FoV is actually far from negligible at these energies. 

3.4. Construction of model lightcurves 

The next step was to construct model lightcurves describing 
the modulation of the radiation of each component described 
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Fig. 3. Model lightcurves of each component for channel 
3 (~27keV) of the first EO. The detector count rate val- 
ues are the effective contribution of the point sources (orange 
stars). The other components are the sky background (red cir- 
cles), the GRXE (blue squares), and the Earth CXB reflec- 
tion (cyan triangles, long-dashed) and the CR-induced emis- 
sion (magenta triangles, short-dashed). For these diffuse com- 
ponents the detector count rate corresponds to an incident radia- 
tion of 40 count s"' sr"' . For the instrumental background (black 
diamonds) we show the relative modulation, Mi„s(f)/Mns, de- 
rived from the SPI/ACS corresponding to a detector count rate 
of 10 counts"'. 



in Sect. 13.21 as induced by the passage of the Earth through the 
IBIS/ISGRI FoV. A complete set of model lightcurves is shown 
inFig.|3] 

3.4.1. Extended components 

For the diffuse components - the CXB, the GRXE and the two 
different Earth emission components - we constructed the model 
lightcurves by generating a series of images of the sky at differ- 
ent times based on the attitude of the spacecraft and the posi- 
tion of the Earth in the FoV (see Sect. 13. II ). For each individual 
component, we considered only its own contribution and the in- 
strumental vignetting effects (see Sect. 13. 3b attenuating the count 
rates on the detector plane. The sum of the pixels in the images 
generated for different times during the Earth occultation defines 
the model lightcurve for a given component. Because the instru- 
mental attenuation is energy dependent, we constructed these 
model lightcurves for each energy bin and for each of the four 
EOs because of the slightly different pointing directions with re- 
spect to the Galactic ridge and Earth positions. As high-energy 
radiation from outside the field of view also contributes to the 
detector counts (see Sect. I3.3l l, the simulated images were de- 
fined on a wide area extending 30°outside of the actual FoV of 
IBIS/ISGRI (|F| < 44.4° and |Z| < 44.6°). 

The normalization of the diffuse components in the input im- 
ages - without vignetting effects - was set to lOcounts"' sr"'. 
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For the Earth emission components, this is the average intensity 
over the Earth disk, while it is the average in the area defined by 
320° < I < 340° and \b\ < 5° for the GRXE. After attenuation by 
the instrumental effects described in Sect. l3.3l the actual detector 
count rate was typically reduced by an order of magnitude (see 
Fig.EJ. 

3.4.2. Point sources 

In addition to the lightcurves constructed for the extended com- 
ponents we also generated one lightcurve in each energy bin for 
the point sources in the FoV. The time modulation is step-like 
in this case due to the abrupt disappearance of a source when it 
gets occultated by the Earth. The count rate used for each source 
detected with a significance of more than 2 cr was derived from 
a bremsstrahlung fit to its observed IBIS/ISGRI spectrum (see 
Sect. |2]l divided by the fraction of time during which the source 
was not occulted. These counts were then assigned to the corre- 
sponding source position in the simulated sky images, and de- 
tector counts were obtained by summing-up the image pixels 
after application of the instrumental vignetting effects. We did 
this at different times during the passage of the Earth through 
the FoV to get model detector lightcurves. As the set of model 
lightcurves for point sources at different energies was based on 
the actual data collected during each EO, they were considered 
as a fix contribution to the detector lightcurves. 

3.4.3. Instrumental background 

The last but the dominant contributor to the observed detec- 
tor lightcurves is the instrumental background. The time vari- 
ability of this component depends on the cosmic particle en- 
vironment and the induced radioactive decay of the space- 
craft's material. The particle environment is well monitored 
by the anti-coincidence shield (ACS) of the sp ectrometer SPI, 
the o ther gamma-ray instrument of INTEGRAL dVedrenne et alj 
12003 ). We found good evidence that the IBIS/ISGRI detector 
lightcurves are indeed following the variations recorded by the 
SPI ACS. To estimate the actual relationship between the count 
rates in the ACS and in the IBIS/ISGRI detector in each of the 
considered spectral bins, we used the extragalactic observations 
of revolution 342 (Her X-1 and XMM LSS). These observa- 
tions, away from bright hard X-ray sources, were taken about 
six months before the EOs and have the particularity of includ- 
ing a solar flare at the start of the revolution, resulting in impor- 
tant correlated variations in the SPI ACS and the ISGRI detector 
counts during the 1 2-hour decay of the flare between INTEGRAL 
Julian dates (IJD = JD - 2 45 1 544.5) of 2040.0 and 2040.5. This 
relationship is characterized by the slope a of a linear fit of the 
ISGRI counts versus the SPI/ ACS counts. This slope is likely to 
change from one observation to the other because the orientation 
of the spacecraft with respect to the solar radiation and particle 
flux will change the effective areas of both the SPI/ ACS and the 
IBIS/ISGRI detectors in a complex manner. However, the energy 
dependence of the slope a(E) for different ISGRI energy bins is 
expected to be rather stable. We used this energy dependence 
of a as an indication of the amount of SPI/ ACS modulation ex- 
pected in the ISGRI detector lightcurves of the EOs. The model 
lightcurve for the variations of the instrumental background was 
thus constructed based on those of the SPI/ACS as: 

MUE, t) = a{E) (CAcs(f) - Cl^) , (2) 
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Fig. 4. IBIS/ISGRI count rate spectra of each model compo- 
nent derived from the detector lightcurves of the INTEGRAL 
EOs. The values are vignetting-corrected count rates per channel 
in units of 10 counts"' sr"\ except for the instrumental back- 
ground (black diamonds) for which they are actual detector 
count rates in counts For the CXB (red circles), the GRXE 
(blue squares) and the total Earth emission (green triangles), the 
dotted lines of the same color show the average spectra obtained 
for the four independent EOs. The relative contributions to the 
Earth emission from the CXB reflection (cyan triangles, long- 
dashed) and the CR scattering in the atmosphere (magenta trian- 
gles, short-dashed) are also shown. 

where Cacs(0 is the SPI/ACS count rate lightcurve measured 
during the EOs. 

3.5. Spectral fitting 

We described above the construction of the model lightcurves 
M,(f) shown in Fig. [3] for each component /, which are 
the SPI/ACS-related variations of the instrumental background 
(ins), the sky background (sky), the GRXE (gal), the Earth's 
albedo (alb), the atmospheric CR-induced emission (atm), and 
the point sources (src) in the FoV. The next step is to adjust these 
model lightcurves to the observed detector lightcurve D(t) in a 
given energy band with a least-square fit. This is done by the 
following linear relation: 

^ins "I" ^ins , 7 "I" ^ins 

'i.:nd ' 

Cgal Mgal(f) + Caib Maib(f) + C^tm M^tm(t) + ^^^(f) , (3) 

where oins, ^ins and the five c, are the seven free parameters of the 
fit, scaling the model lightcurves M,(f) to best match the obser- 
vations D{t). The flins parameter describes the average value of 
the instrumental background, whereas the bi^s and Cins parame- 
ters model its variations. The parameter bi„s allows us to account 
for a linearly increasing or decreasing trend of the instrumen- 
tal background during the observation, centered on a time 7 and 
ending at a time fend- This turned out to have an important effect 
at energies above 60keV (see Sect. l5.4l i. Short-term variations of 
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Fig. 5. Examples of detector lightcurve fits {upper panels) and associated residuals {lower panels) for the first EO at three represen- 
tative energies; ~27keV (channel 3), ~52keV (channel 9) and ~112keV (channel 15), /rom left to right. The reduced value of 
the best-fit curve (red line) is 0.97, 1.08 and 1.06 respectively. The contribution of different components is shown by different colors 
with the same coding as in Figs. [3] and |4] From bottom to top we add to the instrumental background (grey) - modulated by the 
SPI/ACS lightcurve and a possible trend with time - the sky background (red), the GRXE (blue), the albedo (cyan) and atmospheric 
(magenta) Earth emissions, and the fixed contribution from point sources (orange). 



the instrumental background Mins(f) are derived from the simul- 
taneous SPI/ACS lightcurve according to Eq. (O and are scaled 
by the parameter Ci„s. The four other c, parameters are the count 
rates of the various diffuse emission components on the sky, al- 
ready corrected for instrumental vignetting effects. Finally, the 
model lightcurve for the point sources in the FoV, Mstc{t), was 
not scaled as it already represents effective count rates in the de- 
tector. 

By fitting the observed lightcurves D{E, t) in different energy 
bands E, one derives count rate spectra c,(£') for the five com- 
ponents /. When we did this independently for each of the 16 
energy bins, we obtained quite noisy spectra with a divergence 
towards non-plausible values in some channels. This is due to 
significant degeneracy between the various components that we 
discuss in Sect. 15.41 A way to overcome this problem was to 
include a link in the fitting between the values obtained in one 
energy bin and in some others. We did this by adding an addi- 
tional constraint to the;\f^ minimization of the fit so that the fitted 
parameter value c, would not be too far from an expected value 
c™'' according to: 



^^XOogci-iogc;"") 



(4) 



where =X^ l^-oi. is the;if^ divided by the number of degrees 
of freedom (d.o.f.) and ^ is a factor to be chosen to get an ap- 
propriate balance between the quality of the fit and additional 
constraints. We calculated the difference with respect to the ex- 
pected value on a logarithmic scale to avoid a dependence on the 
actual count rate values from one component to the other. 

The choice of the expected values c™'' in this constraint fit 
can of course have strong implications on the results. We there- 
fore took great care to define them without including wrong as- 
sumptions and systematic effects. For the parameter Cins, the ex- 
pected value was set to be the mean value obtained over the 16 
energy bins. This was motivated by our discussion in Sect. l3.4.3l 
where we concluded that this factor can differ from unity, but 
is expected to be rather constant from one energy bin to the 
other. For the four other c, parameters, we defined the expected 



value based on the assumption that the final, unfolded spectrum 
of each component is supposed to be smooth. For each compo- 
nent spectrum, c,(£'), the expected count rate in a given channel, 
c^"'', was set to be the linear interpolation between the values in 
the two adjacent energy bins corrected for the effects of differ- 
ent energy widths of the channels and of the detector respons^B 
For the first and last energy bins, the spectral smoothness was 
similarly constrained, by setting the expected value to the linear 
extrapolation of the two closest channels. We note that we did 
not constrain the spectrum of the instrumental background level 
flins and its increasing or decreasing trend with time, bins, because 
both can change rapidly from one energy bin to the other due to 
the presence of narrow emission lines (see Sect.|2]i. For the five 
other parameters, the introduced interdependence of the values 
obtained in adjacent energy bins typically reduces the number of 
free parameters of the fit by a factor of 2. This was taken into 
account when calculating the d.o.f. of the fit. 

The actual fitting began with a set of input spectra and min- 
imized the modified of Eq. (|4|i, one channel after the other. 
We then reran this process up to four times until the overall - 
computed on the lightcurves in all energy bins - did not signifi- 
cantly improve anymore. This iterative spectral fitting was done 
independently for each of the four Earth-observation datasets, 
and we combined the results to get the final spectra. 

An issue in the fitting process is the choice of the parameter 
^ that defines the strength of the additional constraint on the fit in 
Eq. (01). If ^ is low the spectral smoothness constraint becomes 
important and it then does not leave enough freedom to fit the 
actual data, whereas if ^ is too high, the fit might diverge in some 
energy bins towards unrealistic values. We tested different values 
and chose ^ = 15, which leaves a lot of freedom to the fit, while 



' To take into account the effect of the ISGRI detector response, we 
used XSpec ( Amaud 1996) with the standard ISGRI ARF and redistri- 
bution matrix file (RMF) of OSA 7.0 to derive the relative strength of 
the instrumental modulation of a powerlaw model spectrum from one 
channel to the other. This modulation follows basically that of the ARF, 
with some additional smoothing and distortion towards lower energies 
due to the RMF, with only a slight dependence on the photon index 
taken to be typically of r = 2. 
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limiting strong divergence to only a few energy bins in one of 
the four EOs, namely EO 3 in the 70-100 keV range (see Fig.Hl). 

Another important issue is the choice of the input parame- 
ters, since we experienced that they can have a significant influ- 
ence on the final results. This is due to degeneracy between some 
components that we discuss in Sect. 15.41 It is therefore safe to 
start with values corresponding roughly to the expectations for 
the CXB and also for one of the two Earth emission comp onents. 
We chose the analytical formula of iGruber et al.l d 19991) as the 
basis for the input CXB spectrum. A first guess of the spectral 
values of the other components was obtained by doing a fit with 
the input CXB spectrum fixed and imposing equal contributions 
for the two Earth emission components. We obtained a global 
Earth emi ssion that is much low er at the highest energies than 
derived bv lChurazov et al.l ( l2007h . As the CR-induced emission 
of the Earth is the dominant component at these energies, this 
suggests that its normalization has to be scaled down by a factor 
of ~0.4 (see Fig. [13]). The final set of unperturbated input spec- 
tra was obtained by fitting the lightcurves again, but this time 
with fixed values for both the CXB and the CR-induced Earth 
emission. For the latter, the spectrum was defined by the ana- 
lytical formula proposed bv lSazonov et al ] (120071 Eg. (1)) with 
a normali zation C of 13.2keVcm ~^ s~^ sr~\ i.e. 0.4 times that 
derived bv lChurazov et all (l2007h . 

To smear-out the dependence of the final results on these in- 
put spectra and to estimate the uncertainties, we made a series of 
fits with different input parameter values. We did this by pertur- 
bating each channel value of the initial spectra by a random de- 
viation following a Gaussian distribution with a cr of 30 %. This 
was done independently for all five component spectra, namely 
the CXB, the GRXE, the two Earth emission components and 
the average level of the instrumental background. We performed 
the whole fitting process described above starting from 30 dif- 
ferent sets of perturbated input spectra. For the four EOs, this 
resulted in 120 spectral fits. Instead of taking the average on the 
obtained values, we took the median in each energy bin as a 
robust estimator of the mean. This has the advantage to be in- 
dependent of the use of the actual results or of their logarithm 
and is not influenced by outstanding values. We estimated the 
1-cr (68 % CL) statistical uncertainties on the median taking the 
two results at +34 %/ V4 = +17 % in rank order away from the 
median, where the factor of four stands for the four EOs being 
independent measurements. 

4. Results 

The count rate spectra and uncertainties obtained with the iter- 
ative spectral fitting process described in Sect. [33] are shown in 
Fig. [4] We note that the large scatter from one EO to the other 
is clearly dominating the uncertainties, which suggests that per- 
forming additional EOs in the future will allow us to signifi- 
cantly improve the statistics of the results. We obtained over- 
all average reduced values x]i^a - 1-20, 1.11, 1.15 and 1.13 
for EO 1 to EO 4, respectively. These values only slightly above 
unity show that we get a fair description of all the lightcurves, 
without overinterpreting the data by using too many free param- 
eters. As an example we show the match between the model 
corresponding to the final results for EO 1 and the observed 
lightcurves in three representative channels in Fig. [5] 

The count rate spectra Ci{E) shown in Fig. [4]coiTespond to 
the values before entering the telescope, i.e. they are coiTected 
for the instrumental vignetting effects described in Sect. 13.31 
We could thus directly use them for spectral fitting with XSpec 
to get unfolded spectra in physical units. We did the spectral 




60 80 
Energy ( keV ) 

Fig. 6. Unfolded IBIS/ISGRI spectra of the sky background (red 
circles), the earth emission (green triangles) and the GRXE (blue 
squares) with their best fit (dashed lines) and the more physical 
(dotted lines) spectral models (see Tables [2] and [3]i. The contri- 
bution of the sum of the considered point sources averaged over 
the four EOs is also shown (orange stars). 

Table 3. Spectral fit parameters for the CXB, flie GRXE and flie 
Earth emission. 



Spectra (model") 


n* 






1 2 


N2' 


^red 


CXB (bknpow) 


1.68 


15.2 


28.7 


2.42 




0.51 


CXB (cutoff)'- 


1.95^ 


44. f 


127 






0.85 


GRXE (cutoff) 


0.0-' 


0.43 


8.83 


1.55^ 


0.82 


0.25 


GRXE (bremss)' 




13.9 


14.7 


1.55^ 


0.79 


0.34 


Earth (bknpow) 


0.37 


0.050 


49.4 


1.78 




0.05 



Notes. Model is either a broken powerlaw (bknpow), a cutoff 
powerlaw (cutoff) or bremsstrahlung (bremss), with or without an extra 
powerlaw. Photon index of low- (1) or high-energy (2) powerlaw. 
Normalization at 1 keV in ph cm"- s"' keV"' sr"' . '''' Characteristic 
energy in keV. Either the break energy (bknpow), the cut-off energy 
(cutoff) or the kT energy (bremss). the more physical model. 
fixed parameter value. 



fitting with the standard IBIS/ISGRI ARE and RMF detector 
response files distributed with OS A 7.0. We did not consider 
Galactic hydrogen absorption in the fit because even along the 
galactic plane the hydrogen column density is small enough - 
A^H ~ 2 X 10^2 cm-2 at Galactic coordinates of {l,b) = (330°,0°) - 
to have only a negligible effect. The spectral fit parameters are 
given in Table [3] and the resulting unfolded spectra are shown in 
Fig. [6] with the numerical values given in Table[2] 

The CXB spectrum is best fitted by a broken powerlaw 
model with a break energy at E\j - 28.7 keV and a high-energy 
photon index of T2 - 2.42 ± 0.09, giving a reduced chi-squared 
of xlf.^ - 0.51 for 12 d.o.f. A more physical model for the 
CXB emission - considered as the superimposition of the emis- 
sion of unresolved Seyfert galaxies - is to take a cut-off pow- 
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erlaw model. The degeneracy between spectral slope and cut- 
off energy was solved by fixing the photon index to the value 
of r = 1.95 derived by Beckmann et al. (2009) on average for 
all Seyfert galaxies detected by INTEGRAL. We then obtained a 
good description of the CXB spectrum ;\fj^^^/d.o.f. - 0.85 / 14 with 
a cut-off energy of E - 127 + 20keV, at slight ly higher energy 
than E - 86 keV derived for Seyfert 1 galaxies (iBeckmann et alJ 
I2OO9I) . 

We found that the GRXE spectrum was best fitted by a cut- 
off powerlaw plus a second powerlaw to account for the hard 
tail at energies above ~ 80keV. As the indices of the powerlaws 
are poorly constrained, we fixed their values to Fi = 0.0 for the 
cut-off powerlaw and to r2 = 1 .55 for the hard tail, as derived by 
iBouchet et al.l (l2008h . This gives a very good description of the 
data w ith a;^^^^^ = 0.25 for 13 d.o.f. According to^ Revnivtsev et al] 
(l2006l) . the main population contributing to the low-energy part 
of the GRXE are intermediate polar cataclysmic variables. The 
accretion column onto the magnetic poles of such types of ac- 
creting white dwarfs is emitting optically thin thermal emission. 
The best fit for this more physical bremsstrahlung model is al- 
most undistinguishable from the cut-off powerlaw (see Fig. |6]l 
and gives a typical average temperature of kT - 14.7 ± 1.4. 

The spectrum of the total Earth emission is best fitted by 
a broken powerlaw with a break at £ = 49.4 + 4.9 keV and a 
high-energy photon index of r2 = 1.78 + 0.13. This break en- 
ergy is slightly higher than derive d by th e recent analysis of 
the ^wZ/f/BAT data by lAiello et al.1 (|2008|) . while the obtained 
spectral slope is in remarkable agreement with their result of 
T2 = 1.72 ± 0.08 (90% CL eiTors). We found a different nor- 
malization of the Earth emission spectrum however, which we 
discuss in Sect. 15.31 where we also discuss the separate spectra 
obtained for the albedo and the CR-induced emission. 

The spectrum of the sum of all point sources detected at more 
than 2 o" on average among the 4 EOs is added in Fig.|6]for com- 
parison. The impression that point sources contribute much less 
than the GRXE is misleading. This is related to the arbitrary area 
we chose for the normalization of the GRXE. If we had normal- 
ized it to the area actually covered by the partially coded Fo V of 
IBIS, we would have had a GRXE spectrum scaled down by a 
factor of ~ 4, depending a bit on the EO. This would then lead to 
a higher contribution of the point sources compare d to the GRXE 
in qu alitative agreement with the SPI results bv iBouchet et al.l 
(120081) . 
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Fig. 7. Comparison of the IBIS/ISGRI CXB spectrum ob- 
tained here (red circles) with the previous INTEGRAL results 
of IBIS/ISGRI (black diamon ds), JEM-X (blue squar es) and SPI 
(green triangles) published bv lChurazov et alj (l2007h . 
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Fig. 8. Comparison of the CXB spectrum obtained by 
INTEGRAL - JEM-X measurements (blue squares) from 
iChurazov et al.i ( 120071) and our IBIS/ISGRI results (red circles) 
-with the HEAO-1 spectra and analytical model bv lGruber et al.l 
(Il999h (error bars and dashed line). The data of the A-4 instru- 
ment of HEAO-1 are shown in green with original normalization, 
while we show in orange the spectrum of the A-2 instrument and 
of the model, both increased by 10 % in intensity. 



5. Discussion 

The spectra in the ~ 20-200 keV range presented above will now 
be compared to previously published INTEGRAL results and to 
spectra obtained by other satellites. In the subsections below, we 
discuss this separately for the CXB spectrum, the GRXE and the 
Earth emission. 

5.1. Sky background spectrum 

It is interesting to compare the CXB spectrum obtained by the 
thorough analysis of the IBIS/ISGRI detector lightcurves pre- 
sented here with the IN TEGRAL results previously published by 
IChurazov etaP (l2007h . The comparison is shown in Fig.|7] Our 
approach could significantly increase the useful energy range of 
the IBIS/ISGRI data towards higher energies. The new results 
fall slightly below the previous IBIS/ISGRI spe ctrum, while we 
get a good agreement with the SPI results of IChurazov et al.l 
(l2007i) . except possibly for the first energy bin. 



As shown in Fig. [8] the slightly lower emission we obtain 
now with IBIS/ISGRI is consistent with the HEAO-1 measure- 
ments and its analytical approximation bv Gruber et al.l (Il999h . 
The flux scaling of the HEAO-1 spectrum by -1-10% as sug- 
gested by IChurazov et alJ (l2007h is actually not required any- 
more in the IBIS/ISGRI energy range. The discrepancy appears 
only below 20keV for the INTEGRAL/JEM-X data that indi- 
cate a higher CXB intensity than the HEAO-1 measurements. It 
seems therefore that a simple scaling in intensity of the historic 
HEAO-1 spectrum is not able to consistently adjust the combined 
CXB measurements of INTEGRAL, both below and above the 
turnover 

However, as illustrated in Fig. [H there is some freedom 
within the uncertainties to scale up by ~ 10% the intensity of 
the spectrum of the A-2 instrument of HEAO-1 without changing 
that of the A-4 experiment. This would better match the JEM- 
X measurements and other results by recent X-ray instruments, 
which all suggest a hi gher intensity be low 20 keV than obtained 
by HEAO-l/A-2 (e.g. lGiUi etaIll2007L Fig. 15). The net effect 
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Fig. 9. Comparison of the INTEGRAL IBIS/ ISGRI (red cir- 
cles, this work) and JEM-X (magenta diamonds. IChurazov et alJ 
l2007i) spectra with the other recent CXB measurements by 
Swift and BeppoSAX. The Swift/XRT error box (orange shaded 
area) and the ^wZ/f/BAT results (green triangles) are from 
[Moretti et al. (20( )9) and .Aiello et al. (2008), respectively. The 
best-fit model of iMoretti et al.l (2009) for the combined Swift 
dataset is shown with a black line and grey uncertainty area. The 
original BeppoSAX/PDS measurements of lFrontera et alj ( l2007h 
were scaled by +13% in intensity (blue squares) to correct for 
the difference in Crab normalization with respect to INTEGRAL. 
The analytical model we propose in Eq. Q is shown as a purple 
dashed line. 



would be a broadening of the CXB hump and a slight shift of its 
maximum towards lower energies. The expected qualitative con- 
sequence for an AGN population synthesis of the CXB would 
be a reduction of the contribution of the most highly obscured 
AGN , in particular the Compton-thick ones (e.g. Treister et al.i 
l2009i) . Alternatively, it could also indicate a slightly stronger 
contribution from a population of distant (reds hifted) luminous 
AGN compared to the local population (e.g. iTreister & Urrvl 
I2005h . 

Another consistency check of our results is to compare them 
to the recent Swift and BeppoSAX measurements. This is il- 
lustrated in Fig. |9] where the SwiftK RT and the SwiftfBAT 
specti-a are from More tti et aP (120091) and lAiello et afl (lioOSh . 
respectively. Our data are consistent with the SwiftfBAT re- 
sul ts and the combined XRT-BAT spectral model proposed 
by iMoretti et aP (l2009l) . although they tend to be at a sig- 
nificantly lower inten sity. Our data agree very well with the 
BeppoSAXfPDS data (iFronteraet al.l l2007. Fig. 6 Bottom) pro- 
vided that they are scaled by a factor of 1.13 in intensity to 
account for the diff'erence in the Crab normalization in the 
20-50 keV band between BeppoSAX/PDS ("Frontera et a l.ll2007l 
Foab =9.22x10 ''erg cm^^ g-i) j^^^j INTEGRAL ( Churazo vetall 
12007 . Fcrab = 10.4xl0"''ergcm--s-'). We note that the latter 
estimation for INTEGRAL is fully consistent with the value ob- 
tained with OSA 7.0. The measured fluxes are Fcrab = 10.30 x 
lO^^ergcm"^ s"' and 10.46 x 10"''ergcm"^ s"' for the Crab ob- 
servations of revolutions 365 and 422, respectively. The obtained 
spectrum seems to be also very consiste nt in the peak reg ion 
with the recent CXB synthesis model by Treis ter et al.l(l2009h . It 
thus gives additional evidence for a small Compton-thick AGN 
fraction in the CXB spectrum, close to 9 % instead about 30- 
40 % postulated before (e.g. Treister & Urry 2005). Our data 
cannot constrain a possible hardening of the CXB spectrum 
above lOOkeV, but are consistent with an additional contribu- 
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Fig. 10. Comparison of the obtained GRXE spectrum (blue 
squares) with recent other determinations, all renormalized to 
the central radian of the Milky Way defined by |/| < 30° and 
\b\ < 15°. The previous INTEGRALflBlS data (red cycles) 
and the RXTE/PCA measurements (black diamonds) are from 
Krivonos et al. (2007, F ig. 14). The INTEG RAL/SPl spectrum 
(green triangles) is from lBouchet et a DdaOOSlFig. 9). 



tion of flat-spectrum radio quasars, which have been found to 
dominate the CXB in the MeV range ( Ajello_et al. 2009). 

Based on the considerations above, we can tentatively sug- 
gest a slight adaptation of the analytical description of the CXB 
proposed bv lMoretti et al. (2009, Eq. (4)), as: 

, dNy , 0.109 phcm-^s-'keV-i SI--' 

dE (£/28keV)i-40-H(£/28keV)2-«8 ' ^ 

where the only difiference - but a correction of a typo in the units 
- is a change of the break energy from 29 keV to 28 keV. The 
corresponding spectral shape is at the lower limit of the uncer- 
tainty area of the Swift model as shown in Fig.|9] 

5.2. Galactic ridge emission 

It is not easy to compare results on the GRXE from one pub- 
lication to the other, because the emission is often defined in 
different regions of the Galaxy. As the region covered by our 
observations is away from the Galactic bulge where most deter- 
minations have been made, we have to rescale them to a more 
commonly used area. We choose the central radian of the Milky 
Way defined in Galactic longitude / and latitude b by |/| < 30° 
and \b\ < 15° as the reference area for a comparison of the 
various measurements. As we do have an analytical model of 
the GRXE (see Sect. 13.21 ). it is possible to determine the scal- 
ing factor from any region in the Galaxy to the chosen area. 
The resulting renormalized spectra are compared in Fig.[TOl Our 
measurements had to be scaled by a factor of 0.16 to corre- 
spond to the chosen a rea. The INTEGRAL/IBIS spectrum from 
iKrivonos et"an (l2007l Fig. 14) corresponding to the IBIS FoV 
area centered on the Galactic bulge was multiplied by a calcu- 
lated factor of 2.77. We used the same factor for the RXTE/PCA 
data that have been scaled by Krivonos et al. (2007) to match the 
IB IS measuremen t at 2pkeV. For the INTEGRAL/SPl spectrum 
of iBouchet et al.l (l2008l Fig. 9) that correspond already to the 
area chosen here, we just had to convert the units. 

In general. Fig. [TOl shows a good agreement between the re- 
sults of the various instruments. This is quite remarkable for data 
that were not arbitrarily renormalized, but were rescaled based 
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Fig. 11. Comparison of the obtained GRXE spectrum (blue er- 
ror bars and long-dashed line) with the higher-energy emission 
of the Galactic ridge in the region defined by \l\ < 30° and 
\b\ < 15° as ob served by INTEGRAL/SPl (green triangles) from 
iBouchet et all (|2008, Fig. 9) and by CG/?0/COMPTEL (black 
diamonds) and CGRO/EGRET (red eiror bais). The CGRO data 
are fromjPorter et al. (2008, Fig. 3), which is based on the analy- 
sis bv lStrong et al.i (il999.,2004.) . The green solid Une shows the 
model fitting the SPI observations above lOOkeV (iBouchet et alj 
12008 ^. It is the sum of a powerlaw continuum (green dotted line), 
the positronium continuum (short-dashed line) and a narrow 
electron-positron annihilation line at 5 11 ke V. A scaling by 20 % 
of the high-energy powerlaw is indicated by the blue dotted-line. 



on our very simple double-Lorentzian model of the GRXE. This 
suggests that the model provides a fair description of the overall 
emission of the inner Galaxy in the hard X-ray range. All three 
independent INTEGRAL measurements reveal a minimum at 
about 80 ke V. This wa s only suggested by the 2-cr upper limits of 
iKrivonos et al.l (l2007h . but is confirmed n ow by our Earth occul - 
tation results and the latest SPI results of I Bouchet et all (|2008|) . 
Our IBIS/ISGRI results do, however, suggest that the minimum 
is shallower than previously found. The diffuse GRXE below 
80 ke V is thought to be due to a population of accreting white 
dwarfs too faint to be resolved into discrete sources in the hard 
X-rays (iRevnivtsev et al.ll2006t IKrivonos et"ani2007h . It is only 
at energies of ~ 6-7 keV that the diffuse emis sion could finally 
be re solved using deep Chandra observations ( IRevnivtsev et al.l 
l2009i) . The bremsstrahlung temperature of kT - 14.7 +1.4 
that we derived in Sect. |4] for the accretion column onto the 
pole of the white dwarfs agrees well with the measurements 
of individual intermediate p olar systems detected by Swift/BAT 
dBrunschweiger et al.l2009l) . Based on Table 2 in the latter publi- 
cation, we note that this temperature would correspond to a typ- 
ical whit e dwarf mass of Mw^ ^ 0. 60 ± 0.05 Mq according to the 
model of ISuleimanov et ^ (l2005l) . This estimate is just slightly 
above the expected average mass of white dwarfs in the Galaxy 
(M^d ~ 0.5 Mq) that was found to agree well with the previous 
IBIS/ISGRI results on the GRXE (IKrivonos et al.ll2007l and ref- 
erences therein). 

Above 80keV the GRXE spectrum is likely dominated 
by inverse-Compt on emission from the interstellar medium 
dPorter et al.ll2008l) . We derived an i ntensi ty at the level of the 
2-cr upper limits of Krivon os et al ] (12007), in excellent agree- 
ment with the latest SPI observations (Bouchet et al. 2008). The 
photon index of the high-energy powerlaw derived by these au- 
thors (F = 1.55) also agrees very well with our results, although 
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Fig. 12. Comparison of the obtained Earth emission spec- 
trum (green triangles) with previous determinations by var- 
ious mi ssions. The thin grey line is the INTEGRAL spec- 
trum of IChurazov et al.l (|2007|) as described in Fig. [13] The 
obtained IBIS/ISGRI spectrum lies well between the OSO- 
3 (black diamonds) and the 5w//i'/BA T (r ed circles) measure - 
ments of ISchwartz & PetersonI (11974 and lAiello et alj (|2008|) . 
re spectively. The values of the BeppoSAXfPDS measurements 
of iFrontera et al.l (|2007) were increased by 13 % (blue squares) 
to be consistent with Fig.|9] 



it is too poorly constrained by our data alone to be fitted inde- 
pendently. 

The overall shape of the high-energy spectrum of the GRXE 
including the observations of the Compton gamma-ray observa- 
tory (CGRO) up to lOGeV is shown in Fig.[TT] As our data are 
scaled from a region at 320° < I < 340° lying outside the galac- 
tic bulge where the bulk of the positronium annihilation is emit- 
ted, they should be almost unaffected by the positronium contin- 
uum. This implies that the ~ 20 % higher normalization of the 
high-energy powerlaw suggested by our data should be in t rinsic . 
This would agree well with the discussion of lPorter et al.l (l2008l) 
concerning a possible higher normalization of up to 40 % for this 
powerlaw. 

5.3. Earth emission 

Before discussing the relative contributions of the two Earth 
emission components, we first compare, in Fig. [121 the obtained 
spectrum of the Earth with other determinations. The Earth emis- 
sion is found to be very consistent with the spectra obtained pre- 
viously, although there is a big scatter among the various de- 
terminations. This is at least partially due to the modulation of 
the Earth emission by the solar cycle and a dependence of the 
obser ved flux on the spacec raft altitude and the geomagnetic lat- 
itude dSazonov et al. I I2OO7I) . For instance, the difference in nor- 
malization between the 5w//f/BAT spectrum and our determina- 
tion can be related to INTEGRAL drifting towards an almost po- 
lar orbit, while Swift has a more equatorial orbit. The difference 
by roughly a factor of two depending on the energy is consistent 
with the difference found by the polar-orbiting sate llite 1972- 
076B between the equatorial and the polar regions JTmhof et al.l 
,197 6). We did not include these spectra here for the sake of clar- 
ity, but they would be compatible with the other determinations 
provided that they are corrected for unsubtracted CXB emission 
as shown by Ajello et al. (2008, Figs. 15 and 16). 

A discrepancy we cannot ascribe to a different observation 
epoch or a different viewpoint is the inconsistency of our results 
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Fig. 13. Resulting spectrum of the total Earth emission (green 
triangles, solid line) with separated contributions from the Earth 
reflection of the CXB (cyan triangles, short-dashed line) and 
the CR-induced atmospheric emission (magenta triangles, long- 
dashed line). The thin grey curves are normalized as derived by 
[Churazov et al. (2007), whereas the thick colored lines are nor- 
malized to match our measurements. 



with the spectrum derived from the same INTEGRAL observa- 
tions by Churazo v et al... (2007). We derive a higher Earth emis- 
sion at low energies and a lower intensity at high energies. The 
discrepancy at the highest energies could somehow be due to 
the fact that our results are based on IBIS data and their results 
on SPl, although both instruments are well cross-calibrated. It is 
also possible that the difference comes from our more detailed 
modeUng of several instrumental eff'ects described in Sect. 13.31 
At the lowest energies, the neglection of point source emission 
by IChurazov et al. (2007) is a likely cause of the discrepancy, 
since for a given CXB we need more Earth emission to compen- 
sate for the occultation of point sources. 

To better characterize the difference between the two deter- 
minations of the Earth spectrum, we show in Fig. [T3] the de- 
composition of the overall Earth emission in the two distinct 
components considered in both studies. Those are the reflec- 
tion of the CXB by the Earth - the albedo - and the emis- 
sion induced by CR interactions in the atmosphere. The spec- 
tra of both components can be described by analytical func- 
tio ns fitted to the results of Monte-Carlo simulat ions published 
by Churazov et al.l (l2008h for the albedo, and by ISazonov et al.l 
12007^ for the atmospheric emission. In order to fit the overall 
spectrum of the Earth with these two components, we need to 
increase the albedo component by ~ 40 % and decrease the at- 
mospheric e mission by ~ 60 % com pared to the normalizations 
suggested by IChurazov et al.l (|2007[). We thus get a no rmaliza- 
tion of the atmospheric emission ( Sazonov et al.ll2007l Eq. (1)) 



of 13.2keVcm ^s ' sr ' at the break energy of £ = 44keV, 
instead of 32.9keVcm"^s ' sr ' as derived by IChurazov et al.l 
(I2007h . This difference in normalization is much more impor- 
tant than the few perc ent of likely overestimation mentioned by 
ISazonov et al.l ( l2007l) related to the inclusion of the Compton- 
scattering from particles that would have intersected the surface 
of the Earth. Further speculations on the origin of this discrep- 
ancy are beyond the scope of this paper. 

Concerning the albedo, we note that the increase by a fac- 
tor 1 .4 we derive here is not on the Earth reffection expected for 
the HEAO-1 analytic approximation of Gruber et al. ( 1999), but 
on this spect rum already scaled by a factor of 1 . 1 according to 
the results of IChurazov et al.l (l2007h . As the CXB spectrum ob- 



tained here agrees well with the original HEAO-1 spectrum, our 
results suggest a reffection efficiency of the Earth atmosphere 
~ 50 % higher than ob tained by the Monte-Carlo simulations of 
IChurazov et al.l (l2008h . According to these authors, the shape of 
the input CXB spectrum has only a very limited effect on the re- 
ffected spectrum, especially at energies below ~ 30keV, so there 
must be another reason for the important difference we observe. 
One possibility is related to the delicate modeling of the com- 
position of the atmosphere. Could the presence of clouds have a 
significant effect by increasing the amount of hydrogen atoms in 
the upper atmosphere, resulting in a more " Sun-like" albedo with 
more reflection at the lower energies (see Churazov et alj|2008l 
Fig. 6)? Another possibility would be the presence of another 
Earth emission component emerging at the lowest energies, in 
particular the potentially strong X -ray emission of aurorae (e.g. 
0stgaard et al. 2001). We note that lChurazov et al. (2007) found 
clear evidence for auroral emission in the INTEGRAL/ JEM-X 
data during EO 2 and EO 3. 

Finally, it is fair to mention that the almost perfect agree- 
ment we show in Fig. [13] between the data and the model for 
the two distinct Earth emission components is mainly due to the 
choice of the input spectrum for the atmospheric emission (see 
Sect. I3.5l l. The strong degeneracy between the two Earth emis- 
sion components discussed in Sect. 15. 4! does not allow us to get 
such a good distinction of the two components when starting 
from an arbitrary set of input parameters. We would just get a 
tendency for the CR-induced emission to dominate at the higher 
energies and vice-versa at low energy. 

5.4. Degeneracy issues 

Degeneracy is the reason why it is so difficult to determine the 
CXB spectrum with Earth-occultation data. The basic problem is 
that a spatially uniform Earth emission cannot be distinguished 
from the CXB occultation when the instrumental background 
level is unknown. It is because of this that all previous studies 
using the Earth as occultator had to assume a priori the spectrum 
of the CXB and to a large extend also that of the Earth emis- 
sion (see Churazov et al. 2007; Frontera et al. 2007; Aiello et aTj 
2008). Here, we tried to solve the degeneracy issue by fixing the 
spatial distribution of the Earth emission components rather than 
its spectral shape and by introducing a spectral smoothness con- 
straint as explained in Sect. 13.51 However, because of the noise 
in the data and a possible deviation of the actual instrumental 
background variations compared to those assumed based on the 
SPl/ ACS lightcurve (see Sect. 13.4.31 ) it is not possible to com- 
pletely solve the degeneracy issues. 

To overcome this problem, we had also to incorporate some 
a priori assumptions on the spectral shape of the CXB and of 
one of the two Earth emission components, chosen to be the at- 
mospheric CR-induced emission because of its strong drop at 
low energies that cannot be easily determined otherwise. As ex- 
plained in Sect. 13.51 this is however only used to define the set of 
input spectra that we then perturbate randomly before fitting the 
parameters to the data. Despite this, we still keep a dependence 
on the input parameters in the results. This is very obvious for 
the two Earth emission components that are strongly degener- 
ated (see Sect. 15.3b . 

Another important degeneracy is between the instrumental 
background and the sum of the sky background plus Earth emis- 
sion. A higher instrumental background will imply lower sky 
background and Earth emission, and vice-versa. Actually, the 
data do primarily constrain the difference between the Earth and 
the CXB emissions. This difference is basically the height of a 
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bump in the lightcurve in case the Earth emission dominates or, 
alternatively, the depth of a trough, when the sky background 
occultation dominates (see Fig. |5j- To test the effect of this de- 
generacy on our results, we fitted the data with the same proce- 
dure as explained in Sect. [33] but with an input spectrum for the 
CXB increased by 10 %. This resulted in a CXB spectrum with a 
similar shape, but a higher normalization by about the same fac- 
tor. However, to compensate this, the Earth albedo is then found 
to be higher than derived bv lChurazov et a l. (2007) by a factor 
of ~ 2, instead of ~ 1.4 (see Sect. 15.31 ). Although this cannot 
be completely excluded, the discrepancy on the normalization 
of the albedo is judged to be unrealistically high. We thus favor 
the more conservative results obtained with the original normal- 
ization of the CXB for the input parameters. 

The presence of the GRXE in the observed region of the sky 
can also add some degeneracy, but its location to the side of the 
Fo V was actually rather optimal as its maximal occultation effect 
occuiTed earlier in the lightcurves than for the CXB (see Fig.|3]l. 
We only noted a slight degeneracy between the GRXE and the 
sum of the CXB and the polar-enhanced atmospheric emission. 

We identified an instrumental effect that affects the results 
of both the GRXE and of the Earth emission at energies above 
~ 60keV. By ignoring the possibility of an increasing or de- 
creasing trend of the instrumental background in addition to the 
SPI/ACS modulation (see Eq. (O), we obtained a high-energy 
drop of the Earth emission together with an unrealistically steep 
rise of the GRXE spectrum. This behavior aims at compensat- 
ing a decreasing trend of the instrumental counts in the spectral 
region of the emission lines of W, Pb and Bi (60-80 ke V, see 
Sect. |2| and an increasing trend of the counts at even higher en- 
ergies. This is likely due to radioactive decay at the exit of the 
radiation belts and illustrates the sensitivity of the method on a 
very accurate description of all instrumental effects. 

Finally, we note that an underestimation of the contribution 
of the point sources will tend to increase both the GRXE and the 
CXB intensity because point sources are distributed all over the 
FoV with increased density in the galactic plane. 

6. Conclusion 

We presented the results of an original analysis of the four con- 
secutive Earth-occultation observations by the IBIS/ISGRI in- 
strument aboard the INTEGRAL satellite. Our approach is com- 
pleme ntary to the previous study of these data by Churaz ov et alJ 
( l2007h . because instead of fixing the spectral shape of the CXB 
and fitting its normalization, we attempt to derive the complete 
spectral information from the observed detector lightcurves in 
different energy bins. This requires a deep understanding of the 
instrumental effects and a careful modeling of the spatial dis- 
tribution of the various contributions from the Galaxy, the Earth 
and point sources. Despite inherent degeneracy issues that forced 
us to fit the data with an additional spectral smoothness con- 
straint and adequate input parameters, the approach used here 
results in a coherent set of spectra for the CXB, the GRXE and 
the Earth emission. 

The obtained IBIS/ISGRI results for the CXB are con- 
sistent with the historic HEAO-1 spectrum, without any scal- 
ing in_JntensitY^he_scaling by -1-10% in intensity proposed 
by IChurazov et al.l (l2007h is not incompatible with the actual 
dataset, but is disfavored as it implies a CXB reflection by 
the Earth twice as strong than the one derived by the Monte- 
Carlo simulations of Churazov et al. (2008). The obtained spec- 
trum also agrees well with recent Swift and BeppoSAX deter- 
minations. We propose a slight adaptation of the CXB model 



spectrum suggested by iMoretti et alJ (l2009h . which is based on 
Swift data alone, that implies a reduced fraction of strongly ab- 
sorbed AGN, compared to the HEAO-1 spectrum of Gruber et al] 
(119991) . 

The spectrum of the Earth emission is very well described 
by the contribution of two distinct components: the reflection of 
the CXB that is dominant at lower energies, and the CR-induced 
atmospheric emission. The derived normalizations for these two 
components is however f ound to be very d i fferen t from what was 
suggested by the study of lChurazov et alJ (l2007h . 

With a total observation time of only about a day, these spe- 
cial types of INTEGRAL observations yield a spectrum of the 
GRXE with comparable statistics as obtained by combining all 
available INTEGRAL/SPl observations. This allows us to ob- 
servationally estimate the average mass of white dwarfs in the 
Galaxy. Conducting similar observations in different regions of 
the Galactic plane would be useful to characterize the longitudi- 
nal distribution of the GRXE. 

However, it would be even more important to conduct Earth 
observations away from the Galactic plane to lift any degener- 
acy related to the presence of the GRXE and point sources. This 
would lead to a determination of the CXB with improved statis- 
tics and less systematics and thus fully exploit INTEGRAL'S 
unique capability to observe the entire Earth from a high-altitude 
orbit. 

Acknowledgements. Based on observations with INTEGRAL, an ESA mission 
with instruments and science data centre funded by ESA member states (espe- 
cially the PI countries: Denmai'k, France, Germany, Italy Switzerland, Spain), 
Czech Republic and Poland, and with the participation of Russia and the USA. 
PL has been supported in part by the PoUsh MNiSW grants NN203065933 and 
362/1/N-INTEGRAL/2008/09/0, and the Polish Astroparticle Network 621/E- 
78/BWSN-0068/2008. 



References 

Ajello, M., Costamante, L., Sambruna, R. M., et al. 2009, ApJ, 699, 603 

Ajello, M., Greiner, J., Sato, G., et al. 2008, ApJ, 689, 666 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, 

Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. 

Jacoby & J. Barnes, 17 — b 
Beckmann, V., Soldi, S., Ricci, C, et al. 2009, A&A, 505, 417 
Bouchet, L., Jourdain, E., Roques, J.-P, et al. 2008, ApJ, 679, 1315 
Boughn, S. R, Crittenden, R. G., & Koehrsen, G. P 2002, ApJ, 580, 672 
Brandt, W. N. & Hasinger, G. 2005, ARA&A, 43, 827 

Brunschweiger, J., Greiner, J., Ajello, M., & Osborne, J. 2009, A&A, 496, 121 
Churazov, E., Sazonov, S., Sunyaev, R., & Revnivtsev, M. 2008, MNRAS, 385, 
719 

Churazov, E., Sunyaev, R., Revnivtsev, M., et al. 2007, A&A, 467, 529 
Courvoisier, T. J.-L., Walter, R., Beckmann, V., et al. 2003, A&A, 411, L53 
Frontera, R, Orlandini, M., Landi, R., et al. 2007, ApJ, 666, 86 
Giacconi, R., Gursky, H., Paolini, F. R., & Rossi, B. B. 1962, Physical Review 
Letters, 9, 439 

Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Gruber, D. E., Matteson, J. L., Peterson, L. E., & Jung, G. V. 1999, ApJ, 520, 
124 

Imhof W. L., Nakano, G. H., & Reagan, J. B. 1976, J. Geophys. Res., 81, 2835 
Kinzer, R. L., Jung, G. V., Gruber, D. E., Matteson, J. L., & Peterson, L. E. 1997, 
ApJ, 475, 361 

Krivonos, R., Revnivtsev, M., Churazov, E., et al. 2007, A&A, 463, 957 
Marshall, F E., Boldt, E. A., Holt, S. S., et al. 1980, ApJ, 235, 4 
Moretti, A., Pagani, C, Cusumano, G., et al. 2009, A&A, 493, 501 
0stgaard, N., Stadsnes, J., Bjordal, J., et al. 2001, J. Geophys. Res., 106, 26081 
Paltani, S., Walter, R., McHardy, I. M., et al. 2008, A&A, 485, 707 
Porter, T. A., Moskalenko, I. V., Strong, A. W., Orlando, E., & Bouchet, L. 2008, 
ApJ, 682, 400 

Revnivtsev, M., Molkov, S., & Sazonov, S. 2008, A&A, 483, 425 
Revnivtsev, M., Sazonov, S., Churazov, E., et al. 2009, Nature, 458, 1 142 
Revnivtsev, M., Sazonov, S., Gilfanov, M., Churazov, E., & Sunyaev, R. 2006, 
A&A, 452, 169 

Sazonov, S., Churazov, E., Sunyaev, R., & Revnivtsev, M. 2007, MNRAS, 377, 
1726 



14 



M. Tiirler et al.: INTEGRAL hard X-ray spectra of the cosmic X-ray background and Galactic ridge emission 



Sazonov, S., Krivonos, R., Revnivtsev, M., Churazov, E., & Sunyaev, R. 2008, 

A&A, 482,517 
Schwartz, D. A. & Peterson, L. E. 1974, ApJ, 190, 297 
Smart, D. F. & Shea, M. A. 2005, Advances in Space Research, 36, 2012 
Strong, A. W., Bloemen, H., Diehl, R., Hermsen, W., & Schonfelder, V. 1999, 

Astrophysical Letters Communications, 39, 209 
Strong, A. W., Moskalenko, I. V., & Reimer, O. 2004, ApJ, 613, 962 
Suleimanov, V., Revnivtsev, M., & Ritter, H. 2005, A&A, 435, 191 
Terrier, R., Lebrun, F., Bazzano, A., et al. 2003, A&A, 411, L167 
Treister, E. & Urry, C. M. 2005, ApJ, 630, 115 
Treister, E., Urry, C. M., & Virani, S. 2009, ApJ, 696, 110 
Ubertini, P., Lebrun, F, Di Cocco, G., et al. 2003, A&A, 411, L131 
Ueda, Y., Akiyama, M., Ohta, K., & Miyaji, T. 2003, ApJ, 598, 886 
Vedrenne, G., Roques, J.-P, Schonfelder, V., et al. 2003, A&A, 411, L63 
Winkler, C, Courvoisier, T. J.-L., Di Cocco, G., et al. 2003, A&A, 411, LI 
Zdziarski, A. A. 1996, MNRAS, 281, L9+ 



